Gluon chain formation in presence of static charges 
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We consider the origins of the gluon chain model. The model serves as a realization of the dynamics of the 
chromoelectric flux between static quark-antiquark sources. The derivation is based on the large-A'c limit of the 
Coulomb gauge Hamiltonian in the presence of a background field introduced to model magnetic confinement. 



I. INTRODUCTION 

The gluon chain model of Greensite and Thorn iden- 
tifies the chromoelectric flux tube that exists between static 
quark charges with a string of quasi particles, constituent glu- 
ons. Through lattice simulations and phenomenological anal- 
yses it is well established that the instantaneous. Coulomb po- 
tential between static charges is confining ||5]-(8l. Even though 
it does not correspond to a physical observable, the static po- 
tential does provide physical insight into the possible origins 
of the confinement mechanism as illustrated by the Gribov- 
Zwanzinger model |9, 10] and other, e.g. variational mod- 
els [1 1- 15,1 . Lattice simulations indicate that the correspond- 
ing string tension is larger (by a factor of 2 to 3) as compared 
to the string tension extracted from time-dependent large Wil- 
son loops. This is consistent with expectations of variational 
analysis. At fixed quark-antiquark separation the Coulomb 
potential corresponds to the energy of a quark-antiquark pair 
in a vacuum state that is unmodified by the presence of the 
pair while the energy extracted from the Wilson loop corre- 
sponds to the energy of the exact QCD eigenstate in which 
the quark-antiquark {QQ) pair polarizes the vacuum lfT6l . The 
gluon chain model is a particular realization of the latter, i.e., 
the exact pair state. Confinement originates from the conden- 
sation of chromomagnetic charges llT7l - l20l . Formation of the 
gluon chain should therefore also provide insights into the in- 
terplay between constituent gluons and magnetic domains in 
the vacuum. 

In the Hamiltonian formulation the true QQ state is gen- 
erated by the evolution operator limp^oeexp(— )3//) from the 
unperturbed vacuum. This is because in a physical gauge the 
Hamiltonian H contains all gluon interactions which also cou- 
ple to the classical, external quark-antiquark color source. In 
this paper we investigate if/how the gluon chain emerges from 
the evolution operator. We follow a canonical formulation of 
QCD in the Coulomb gauge since it contains only physical 
degrees of freedom, and these can be directly related to quasi 
particles. The gluon field is decomposed into normal modes 
representing particle excitations, and a physical state is repre- 
sented as a superposition of multi-gluon states. Furthermore 
the normal mode expansion can be performed with respect to a 
non-vanishing classical background. Such a background is in- 
troduced to (phenomenologically) parametrize topologically 



disconnected sectors of the vacuum. In terms of the path inte- 
gral representation these sectors correspond to large field con- 
figurations, i.e., field domains that cannot be smoothly con- 
nected to the null field configuration |21 1. 

The paper is organized as follows. In the next section we 
review the structure of the Hamiltonian, introduce the particle 
basis, and discuss the role of the individual interaction terms 
in the formation of the chain. In Sec. Ill we propose a sim- 



plified computational scheme for studying the formation of 
the chain state and discuss numerical results. A summary and 



outlook are given in Section IV 



II. QCD HAMILTONIAN AND GLUONS 

In the Coulomb gauge |22| the gluon field is described 
by the vector potential, A"(x) that, for each color compo- 
nent a — 1---Nq— 1, satisfies the trans versality condition, 
V • A" = 0. In the Schrodinger representation the conju- 
gate momenta, which are proportional to the electric field, 
are given by 11" (x) = — /5/5A"(x). The temporal compo- 
nent of the gluon field is eliminated using Gauss's law. This 
leads to an instantaneous interaction between color charges. 
The total color charge density has two components, p(x,fl) — 
Pi,{x,a) + Pg{x,a), coiTesponding to gluons and quarks, re- 
spectively. In the following we ignore dynamical quarks, and 
the only quark charge we consider is that of a static quark- 
antiquark pair placed along the z-axis a distance R apart. The 
corresponding density is therefore given by 

p,(x,fl) = Q]{x)TrjQj{x) ~ Q]{x)Tf,Qj{x). (1) 

Here Q] {x){Qi{x)) represents an operator that creates (anni- 
hilates) a quark at x in a state with color i — 1 • • -Nc, and T" 
are the SU{Nc) color matrices in the fundamental represen- 
tation. We suppress the (iiTelevant) spin indices. Similarly 
QI{x) and Qi{x) are the creation and annihilation operators 
for antiquarks. A state with a static QQ pair is created by the 
operator Qj {zR/2)Qj{—zR/2). The gluon charge density is 
given by 



p,(x,fl)=/„,,A'^(x).n'-(x), 



(2) 
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and the Hamiltonian takes the form 

H = Hfc +Hb+Hc 
where the kinetic plus magnetic terms are given by 
1 



(3) 



Hk+Hb = - I dxi J-'[A]nj-'[A]n + B'), (4) 



and 



He = - J dxdyJ-'[A]pix,a)J[A]K"\x,y,[A])piy,a) 

(5) 

represents the instantaneous Coulomb interaction between 
color charges. Here ^[A] — Det(--D • V) is the Faddeev- 
Popov determinant, D — D,,^ = ^Sab + gfacb^'^ is the covari- 
ant derivative, and B = B" = V x A" +gfabc^^ x A" /2 is the 
magnetic field. The non-abelian Coulomb kernel is formally 
given by 



/:(x,y,[A]) = (D.V)-i(-g2y2)(D.V) 



(6) 



The above describes the Hamiltonian in the Schrodinger rep- 
resentation. The particle basis representation is obtained via 
a canonical transformation from A, 11 to a set of operators 
a^(k,A,fl), a(k, A,fl) representing creation and annihilation 
of gluons with three-momentum k (A; = |k|, [dk] dk/(2;r)^), 
heUcity A, and color a 



A"{x 



[dk] 



1 



e(k,A)a(k,A,fly'^-'' + /i.c] 

n"(x) = - [[dk]J^[Y, e(k,A)a(k,A,fl)e*-''-/;.c], 
I J V 2 



(7) 



Particle operators satisfy ladder algebra and generate a Fock 
space labeled by the number of gluons, n,, occupying a state 
of a given momentum, helicity and color, / = (k,A,fl) 

|«,,«2, •••«,■•••) - {air • • • («/)"' • • • |0)- (8) 

The state with no gluons, |0) = |0,0, • • • ) is annihilated by all 
annihilation operators a, . 



A. The vacuum state 

In the absence of quark sources, after normal-ordering the 
gluon operators, the Hamiltonian 



H={0\H\0)+:H: 



(9) 



contains an infinite number of terms that connect states with 
any numbers of gluons IITTIl . The ground state, can there- 
fore be formally written as 



in) = 



LI 

"1 f^i 



(10) 



The non-uniqueness associated with the definition of a gluon 
state, and thus the Hamiltonian in Eq. (|9|, arises from the ar- 
bitrariness in the choice of the function (o{k) in Eq. (|7]l. For 
example the choice Co{k) = k corresponds to a basis of non- 
interacting particles which diagonalizes the free Hamiltonian 
(i.e., for g = 0). Other proposals, based on the variational 
principle, have been analyzed in fTP-TSl. These studies con- 
sidered an optimal choice for the basis of states obtained with 
(0{k) that approaches the free particle limit for large k and 
is large and possibly divergent in the infrared (IR), i.e., for 
— > 0. This is because an IR enhanced Co{k) suppresses con- 
tributions to vacuum expectation values (vev) from fields near 
the Gribov horizon 1 23 1 and removes the Landau pole from the 
Coulomb kernel (cf. Eq. (j6]l). With such an optimal choice the 
vacuum in Eq. (10 1 is approximated by the state with a van- 
ishing number of gluonr^ i.e., ~ |0), and the ground state 
energy is therefore given by the first term in Eq. (|9]). 



B. The variational QQ state 

We next consider a state containing the QQ pair. A vari- 
ational state, |/?), which does not take into account the back 
reaction of quarks on the vacuum can be defined as (in the 
volume y) 



^Gj(^£)ej(-f£)|o), 



(11) 



and it is normalized, (R\f() = 1. Even if |0) was the exact 
ground state, this state would only be an approximation to 
the exact QCD eigenstate containing the QQ pair. This is 
because with Pq the term in He proportional to p^ x Pg 
does not conserve the gluon number. The expectation value 
of the Hamiltonian in the variational QQ state defines the 
Coulomb potential, Vc{R), which is proportional to the expec- 
tation value of the Coulomb kernel in the variational vacuum. 



Vc{R)5ab^^{0\K"\R)\0) 



(12) 



Here K"\R) is given by Eq. (|6| evaluated at the positions of 
the quark and the antiquark. The vacuum expectation value 
may be computed by expanding the covariant derivatives in 
powers of A (cf. Eq. (j5]l) and noticing that in the variational 
vacuum. 



(0|A«(x)A^(0)|0)=5, 



'ab 



dk g7-(k) 
(27r)3 2(o{k) 



(13) 



where 5r(k) = I;L=±e,(k,A)e*(k,A) = 5ij-kikj/k^. The 
behavior of Vc{R) at large-7? is correlated with the IR behavior 
of co{k). While early variational studies indicated that with a 



' More accurate approximations, whicli take into account residual correla- 
tions among the "optimal gluons." can be constructed using the standard 



many-body techniques of cluster expansion I24II25I . 
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proper choice of (o{k) it would be possible to obtain a confin- 
ing potential, more detailed analyses showed that all solutions 
are massive, i.e., when transformed to momentum space Vc{k) 
is always finite in the limit A: — > 0, a.k.a. non-confining flSl. 
We now believe this is consistent with lattice results. As 
shown in Q the large-^ strength of the Coulomb potential 
originates from magnetic charges in the vacuum. These are 



absent in the variational model calculation of Eq. ( 12 1 that 
is driven by fields in the neighborhood of the A config- 
uration. This is because magnetic charges are topologically 
disconnected from the first Gribov region where the expan- 
sion applies. Thus it is likely that the string tension, Oc, of the 
variational model Vc{r) ^ Ocr should at most only be a frac- 
tion of the Coulomb string tension and, more likely, Vc of the 
variational model ought not to be confining. In the following 
we further explore these scenarios. 

It is straightforward to show that the expectation value of 
the Hamiltonian in the variational QQ state is given in terms 
ofVcby 

{R\H\R) - (0|//|0) = CfV,{R)-CfVM (14) 

where the last term arises from self-energies of the two static 
quarks {Cp = {N^ - \ )/2Nc is the SU{Nc) color Casimir in 
the fundamental representation). As already mentioned above 
the Coulomb term. He, involves coupling between quark and 
gluon charges. It seems reasonable to expect that this interac- 
tion might be responsible for generating the gluon chain. In 
the particle basis the gluon charge density is given by 

Pg(x,fl) = I^P/'(x,fl)a/a,+J^(p5(x,fl)a/a]+/!.c.). (15) 

The first term is diagonal in the particle basis and because |0) 
contains no gluons it vanishes when applied to the QQ state 
defined by Eq. ( [TT| . The second term, however, changes the 
number of gluons by two and thus could be generating the 
chain. We will return to this possibility below. There are other, 
more complicated interactions involving the quark charge and 
gluon operators that change the number of gluons. They orig- 
inate from the A-dependence of the Coulomb kernel. In the 
particle basis, the Coulomb kernel can be written as 



K'"'{R) = ~V,{R)+ : K'"'{R) : 
where the normal-ordered part is given by 



(16) 



{«},{m} 



,,„,... (a/)'"(«lr---arar---. 

(17) 

Here A'|„} |,„} are the matrix elements of the full kernel eval- 
uated between states containing {«} and {m} gluons, re- 
spectively. Thus, when multiplied by the normal-ordered 
Coulomb kernel mixes the variational QQ state with states 
containing arbitrary numbers of gluons. As shown in ll26l . 
however, in the large-7? limit the matrix elements ^Tj^j 
for {OT},{n} ^ are expected to be smaller than those for 
{m} = {«} — 0. Therefore we expect that at large-/? the dom- 
inant interaction between quark sources and dynamical glu- 
ons originates from the off -diagonal gluon charge density {c.f. 




FIG. 1. Interaction between quark charge (upper line) and the off- 



diagonal gluon charge pfj which creates two gluons 



The dashed 

line represents the Coulomb potential given by the v.e.v of the fully 
dressed Coulomb kernel. 



Eq. ( [T5| )) coupled to the quark charge via Vc, and is given by 



flfxflfyp;' (x)y, (|x - y I ) ^(p5 (y , fl) a/ a] +h.c), (18) 



and shown in Fig. ([T]). In Eq. ( 18 1 the gluon charge density 
creates (annihilates) two constituent gluons in a color anti- 
symmetric state. Thus the combined spin and spacial wave 
function of the gluon pair also has to be antisymmetric. How- 
ever, since Pg is a scalar under rotations, the matrix element, 
pfj is symmetric in spin and relative momentum. Thus the 
above candidate operator for the gluon chain actually vanishes 
identically. 

The variational basis based on the mode expansion in 
Eq. (|7]i seems incompatible with the gluon chain picture. 
There is further evidence that a model in which the vacuum 
is described solely in terms of fluctuations around the A — 
configuration, as implied by Eq. (j7]l, is inadequate. If Vc is 
confining then the expectation value of // in a single gluon 
state is infinite |,27J at all temperatures, and the model fails 
to predict the deconfinement phase transition ||28l . It is well 
established that confinement is related to the presence of mag- 
netic domains in the vacuum, and these are absent in the varia- 
tional vacuum state. One would expect that the magnetic term 
should play an important role in confinement since even 
the classical Yang-Mills field equations have monopole solu- 
tions Il29l . 

In the presence of QCD instantons (a.k.a. monopoles) quan- 
tization has to be performed in each topological sector In our 
phenomenological approach we approximate this by general- 
izing the mode expansion of Eq. (|7]i to describe field fluctua- 
tions, Af, with respect to a classical background field, Ag. 



A"(x)^A«^(x)+A^(x). 



(19) 



This classical field mocks the nontrivial topological vacuum 
and will be specified later Thus Eq. dTli now applies to 
A^ = A — Ab and Ily = 11. Since Eq. Jl9]l is a canonical 
transformation the Hamiltonian can be obtained by substitu- 
tion. Thus in the background field, at large-R, the dominant 
contribution to the Coulomb interaction between quark and 
gluon charges is given by 



He 



(20) 
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Here i/^^y is the interaction between quark charges mediated 
by the Coulomb potential, 



qq 



t/xt/yp^(x,fl)K(|x-y|)p^(y,fl), (21) 



H^g is the quark-gluon charge density interaction diagonal 
with respect to the gluon number, 



H: 



dxdypJx,a)V,{\x-y\)p^{y,a) (22) 



with = Y^ipj {x,a)al ai, and H^^ is the normal-ordered, 
diagonal interaction between gluon charge densities 

Hg = -:J fl'xfl'ypf (x,fl)y,(|x-y|)pf (y,«) : . (23) 

Finally the two terms proportional to Ag, //^ and i/^, are 
given by 



H„h + H: 



M _ 



J dxdyp^(x,a)V,{\x-y\)p^{y,a), 



with 



p„^(x,fl)=/„,,A*(x)n^(x) 



(24) 



(25) 



and describe the interaction of physical gluons with the back- 
ground field and the gluon pair creation in the presence of 
the background, respectively. Physical states should be color 
neutral, thus creation or annihilation of a single gluon can be 
neglected. In the presence of the background, the expectation 
value of the charge operator 



Qi=UJ dxA's{x)n'{x) 



(26) 



in physical states vanishes. However, in a simple classical 
model for the distribution of background fields, as described 
in Appendix |A] quantum charge fluctuations do not vanish, 
i.e., Q"Q" 7^ even for color singlet states. We thus modify 
the right hand side of Eq. (|24]l in such a way that these fluctu- 
ations do not contribute to the energy, yielding 

V j dxdyp^{x,a)Vci\x-y\)p^{y,a) 

+ Vc{0)Q"bQ"b. (27) 



H„h + H, 



After normal-ordering, the term in Eq. (24i proportional to 
a^a defines and the term proportional to a^a^ +li.c. 
gives The difference between the gluon density-density 



interaction and the normal ordered Hamiltonian of Eq. (23 i is 
proportional to either a^a or aa +h.c. These, together with 
the kinetic and magnetic terms combine to ifTTl i) renormalize 
CO via a gap equation which eliminates terms proportional to 
aa + h.c, and ii) modify the single gluon energy. Thus the 
final Hamiltonian can be expressed in the form 

= Y,Eia}a,+H,,+H^g+H^g +H^^+Hf^ (28) 

where £, = E (k) is the single gluon energy in the presence of 
the background field. The action of these operators on gluon 
chain states is shown in Figs. (j2]), (|3]l, Q, and (jSj). 



FIG. 2. Illustration of the matrix element of the one-body operator 
Hamiltonian, Eia'a,, in the gluon chain state with A' = 1. 
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FIG. 3. Part of the leading contribution at large-A'c from //^^ which 
corresponds to quark or antiquark self energy, shown here for the 
A' = 1 gluon chain state. 
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FIG. 4. Quark-gluon, /f^, and gluon-gluon interaction, H^g, matrix 
elements, for A' = 1 and N = 2 chain states, respectively 
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FIG. 5. Diagonal and off-diagonal //^ interactions between 
quasi-gluons and the background field (shown as blobs), for N = 3 
gluon chain states. 



C. The basis for the gluon chain 

We define the chain in a large-A^ limit by a model in which 
the gluon chain state is a superposition of multi-gluon states. 



\QQ.R) 



N 



(29) 



with each state in the sum describing a product of single 
gluons ordered in color and space along a straight light be- 
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tween the quark-antiquark sources, 



nn- I '^^N I dxN^i--- / dxi 
Ty/Nc J-R/2 J-R/2 J-R/2 



R 



Qli^z [G\x^)---G\x,)]uQ]{~-z \0). 



R 



(30) 



In the large-/? Umit the longitudinal, i.e., along the QQ axis, 
and perpendicular motions of gluons factorize. The spacial 
distribution of gluons in the plane perpendicular to the QQ 
axis is given by a single-particle wave function, yf{k,X) — 
e*'(k, A)v/'(kx), which defines the gluon operators, G, in the 
chain (x — {Q±,x)) 

4W =L / [^k]a+(k,A,fl)7;.Xk,A)e-''^-^. (31) 

The normalization constant is obtained from (A^|A^) = 
Zjj{CFlR)^ /T{N + 1) = 1 where / is the normalization inte- 
gral for the spacial wave function, Xj/, ([dk±] = d^kj_/ (2;r)^) 

/ = (wlw) = /[«'kx]v/'(kx)4'(k±)r^'(kx). (32) 



In the large-A^c limit, computation of the leading contributions 
to the matrix elements of the effective Hamiltonian of Eq. (28 i 
in the basis of the gluon chain states, Eq. ( [30| , is straightfor- 
ward. The details and numerical results are presented in the 
next section. 



III. FORMATION OF THE GLUON CHAIN AT LARGE 
QQ SEPARATION 

As discussed in Sec. lIIBI one could consider two models for 
Vc{R). In what we refer to as model-/ Vc{r) will be linearly 
confining and of the form 



y/(r) = CJ,r + y,(0), 
and in model-// Vc is asymptotically flat, 

\imV^.'{r)^Vc{oo) <« 



(33) 



(34) 



We concentrate on the interactions induced by the effective 
Hamiltonian in the limit of large quark-antiquark separation. 



A. Matrix elements of the effective Hamiltonian in the chain 
basis space 



The one body term, Hg, in Eq. ( 28 1 acts independently on 



individual gluons in the chain created by the operators G' (cf. 
Eq. (|3T])). Using 



we find 

(A^|//,|A^)=z2(Cf(v/|v/» 

xCf{w\E\w) 
,{W\E\W) 



N 
i=l 



R/2 
-R/2 



dx, 



dxi 



R/2 



-N- 



{v\v) 



dxi^ 1 

R/2 

-:N-^\e~ 



R/2 



dxi 



VciO)] 



(36) 



where 

{W\E\W) 



j[dki 



]£ ( |kx I ) V^' (kx ) 5'^ (kx ) (k± ) • (37) 



and to define e we subtracted from the single gluon energy 
a constant proportional to the negative of the potential at the 
origin. In color singlet states the total energy of the system 
should be invariant under a constant shift ll30l[3Tl . which we 
now demonstrate. The single gluon energy, E{k), contains self 
energies. In the variational approximation the component of 
the self energy due to the Coulomb interaction is given by ifTTI 



Nc 



[fl'qlK.Ck-q) 



1 + k-q to(^) 
2 0}{q) 



(38) 



where Vc is the Fourier transform of the Coulomb potential. 
For a linearly rising, confining potential, e.g., model-/, the 
low momentum singularity of V{k) is not integrable and the 
resulting infinite self energy can be interpreted as a manifes- 
tation of confinement of color charges. A finite self energy is 
obtained by subtracting the IR singularity which leads to 



with finite and given by 

that follows from 



Nr 



y,(0) 



1 + k-q a)(^) 
2 (o{q) 



-1 



y,(0) = / [dq]VM). 



(39) 



(40) 



(41) 



Even though for a confining potential the Fourier transform 
is defined modulo a constant, it is expected that when all, 
self and mutual, interactions between color charges are ac- 
counted for the dependence on ^ (0) disappears from color 
singlet matrix elements. This will also be the case for the ma- 
trix elements of the effective Hamiltonian in the chain basis 
considered here. In anticipation of this result, in Eq. ( [36| we 
defined an IR finite single particle energy e{k) by subtract- 
ing the Coulomb self energy equal to — A^fV^ (0)/2. Thus, in 
the last line of Eq. (|36]), e is finite, and the IR singularity of 
the confining Coulomb potential is explicit in the term pro- 
portional to Vt(0). In the case of model-// with non-confining 
interactions, self-energies are IR finite but we can perform the 
subtractions nevertheless. 

In the absence of chained gluons, N = 0, the interaction 
between quark densities produces the Coulomb potential be- 



[G{x),G'^ {y)]ij — CF{yf\w)5{x — y) = Cf/5(x — y), (35) tween quark charges (cf. Eq. ^T4\). With gluons separating 
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FIG. 6. A non-planar diagram, which we do not take into account, 
representing direct Coulomb interaction between the quark and the 
antiquark in the presence of 2 gluons. 



the quark from the antiquark, the direct interaction between 
quark charges is nonplanar (c/. Fig. (|6]l) and suppressed by 
a power of Nc compared to successive Coulomb interactions 
between the quark and the nearest gluon or the interaction be- 
tween any two nearest-neighbor gluons in the chain. To lead- 
ing order in the H^iq contribution thus reduces to the quark 
self energies. 



{N\Haq\N) ^ ^CfVM 



Nc 



(42) 



The quark gluon interaction, to leading order in Nc, couples 
the quark (or the antiquark) to the nearest gluon in the chain. 
For example, for the antiquark-gluon interaction we find 



{N\H^^\N)^zUCF{w\w)f-' 

rR/2 



X / dxN ■ 

l-R/2 



-R/2 



dx,CF{gQ\H,„{xi)\gQ) 



(43) 





:VA'(kx)[5r(kx)5r(qx)]'V^'(q±)- 



In the limit R ^ °° where xi/R = 0{l) this reduces to 



(44) 



(^e|//,,(xi)|ge) = ^y,(^|^-xi|^(V/|v/). (45) 

Taking into account both quark and antiquark contributions, 
for the H^„ matrix element we obtain 



{nK\N) ^NNc f'\z (z+ 1)" V. [r {z- 1)) . 

(46) 



For the linearly rising potential of model-/ Eq. (46 1 yields 



N+\ 



-NcVc{Q). 



(47) 



while in the case of model-// we find 



{N\H^ANf=NcV,H. 



(48) 



The interaction between two nearby gluons in the chain given 
by H^g is also straightforward to compute, and passing di- 
rectly to the R^ oa limit we find, 

N-\ nRjl 

{N\H^^\N}=zj,{CF{ilfmf-'L dx^--- 



f-W n+1 fNc 

X / dxi+i / dxiCf — 

J-R/2 J-R/2 V 2 

X / dxi^i ■ I dx\ 

J-R/2 J-R/2 



Vc{xi+i -Xi){\if\\ir)^ 



Nl /-i/z 

Cf J-l/2 7-1/2 



.»,f)%,.,.,.-.„<l±^. 

(49) 



For the linear potential of Eq. ( 33 1, to leading order in Nc this 
yields 

{N\Hf^\Ny = ^^acR + ^iN-l)VM (50) 
and for the asymptotically constant potential 

{N\H^^\N)" = ^{N-l)V,H. (51) 



Since all terms in the effective Hamiltonian (including the self 
energies) are 0{g^), and limNi^.^oo Ncg^ = 0{l), at large Nc 
all matrix elements are finite when expressed in terms of e = 
Nce/3, Gc = Nc(7c/3 for model-/ and K.(°o) = NcVc{°°)/3, 
Vc{Q) = NcVc{0) /3 for model-//, respectively. 

Adding all diagonal contributions of the effective Hamilto- 
nian matrix that are independent of the background field, we 
thus find. 



and 



{N\H\Ny ^^Ne+^acR 



{N\H\N)"^^Ne+^{N+l){V,H 
3 

= -Nm„ + c 



(52) 



-VM) 



(53) 



for model-/ and model-//, respectively. For N ~0 this agrees 
with Eq. ( [I4) i, while, for N> 1, eigenstates of Eq. ( 52 1 or ( 53 1 
represent a tower of chain states with energies proportional to 
the number of gluons in the chain. Clearly the lowest energy 
state of the diagonal part of the Hamiltonian is the variational 
QQ state, with N = gluons. The genuine chain contribution 
to the lowest energy state must therefore originate from the 
terms in the Hamiltonian which couple the constituent gluons 
with the background field, as expected. The interaction of 
physical gluons with the background is given by 
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N-l pR/2 rXi+2 

{N\H^,\N) = zl{CF{xif\xi/) f dxN--- dx,+i 

* ~{ J-R/2 J-R/2 



which gives 

{N-2\H^,\Ny={N\H''\N-2)' = 



3^N{N-i: 



2N{N+l) 



X dyi dxiFB{\yi-Xi\) dxi^i--- dxi 

J-R/2 J-R/2 J-R/2 J-R/2 



(62) 



and 



(54) 



{N-2\H%\N)" = {N\HfAN-2)" 



where 



IrRj^iVcH-Vcio))- (63) 



FBi\y-A) 



Nc 



Y{\y-z\){Vc{\y-z\)-Vc{Q)), (55) for the two models, respectively. 



and 



7 = y Jxxfl'yx[a'kx][«'q±]V«(k±)«(q±)e-''»^''^+*^^ 
[r(q±)5r(q)]'Gy(xx-yx,3'-z)[r(kx)5r(k)]^- 



ivlw) 



(56) 



The correlation function Gc is obtained from the density of the 
vacuum fields 



G'^'(xx-yx,x-y) 



(A-(xx,x)Af(yx,3')> 



1 



(57) 



Here the expectation value is taken with respect to the dis- 
tribution of sources of the background field. These might ef- 
fectively describe monopole-antimonopole pairs in 3D, vortex 
surfaces in 4D, merons, etc.. A simple model is considered in 
the Appendix. Since it is these background fields that are re- 
sponsible for confinement in the first place, i.e. generation 
of the Coulomb potential V<; (^), we assume that the density 
of the underlying magnetic sources is approximately uniform 
over the quark-antiquark separation. So for |x — y| <R we 
expect in general 



G'-''(xx-yx,x-3')^G'''(xx-yx) 



(58) 



and 7 in Eq. (55 1 reduces to a constant of 0{Aqcd), i-e. it is 
independent of the longitudinal distribution of gluons along 
the chain. 

For model-/ evaluation of the integrals in Eq. ( [56| then 
gives 

\ ^ gbi / t Y{N + 3)' {N+l){N + 2)' ^ 

(59) 

while for model-// we find 



N-l 
N+l 



(K.H-K.(o)). 



(60) 



Finally we consider the components of the interaction be- 
tween physical gluons and the background that changes the 



gluon number. From Eq. (24 1 we find (for N >2) 



{N-2\Hf,\N) = {N\Hf^\N-2) = ZnZ^_2{Cf{wW)) 



N-\ 



N-l ,.R/2 



dxiFsixi+i -Xi)--- dxi 

-R/2 J-R/2 



(61) 



Collecting all the terms, Eqs. (52 1,(59 1,(62 1 for model-/ and 



Eqs. ( 53 1,( 60 1,( 63 1 for model-//, we find the following expres- 
sion for the matrix elements of the Hamiltonian in the gluon 
chain basis for large-A^, 



{N'\H\N)' = ^NedN'N+lR^c (l +'-^r) 

3 R 3 R 

+ 2y^cj^RSN'.N-2 + 2y^cj^RSN'-2.N (64) 



(N^lHlNf - -NmgdN,N+-ryR{VcH-V,{0))SN'N 

+ IrRiVeH - Vcm{5N',N-2 + Sn'-2,n) 

(65) 

Here r is the ratio of the diagonal to off-diagonal matrix ele- 
ments in the limit of large-/?. The specific value r = 2 follows 
from the fact that in the two models both terms originate from 



the same interaction cf. Eq. ( 27 1. Below, while presenting 
numerical result, we will also discuss the dependence of the 
lowest eigenvalues on this ratio. 



B. Numerical Results 

Before analyzing the spectra of the effective chain model 
Hamiltonians we consider the large-/? limit of the matrix 



{N'\H\N)^R{5n',n-2 + Sn'-2.n)- 



(66) 



It is straightforward to show that the ground state energy of H 
for large-/? is —2R. For the Hamiltonian of model-/ this im- 
plies that if the kinetic term (proportional to A^e) was ignored, 
the lowest eigenvalue of H' for large-/? would behave as 



3 _ / R 

2^^<^+'"M 



(67) 

At r = 2 the quadratic term vanishes and the lowest chain 
state energy is expiated to grow linearly with /?. At large-/? 
if r >> 2 then the lowest eigenvalue is dominated by the di- 
agonal term. In this case the expectation value of A^, 



(68) 
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FIG. 7. Ground state energy (solid line) of the chain Hamiltonian 
in model-/. A power law fit yields = 0.984(/?GeV)°-^^^GeV. The 
dashed line gives the energy of the bare state using for the string 
tension dc = 0.1 GeV^. 

where y/o is the wave function of the lowest energy chain state, 
can be determined by minimizing the diagonal part with re- 
spect to A^. This gives 




and the ground state energy approaches 

El, = 3R^rYace+^Rac. (70) 

Thus for r > 2 the energy of the chain is higher than the en- 
ergy of the bare state, |0). If r < 2 the off-diagonal term dom- 
inates and the ground state energy becomes negative and pro- 
portional to —R^ while the average number of gluons in the 
chain (A^) —0{l). However, when the kinetic term is included 
in the critical case r = 2 the lowest energy of the chain state 
no longer increases linearly with R. After numerical diagonal- 
ization we find 

for a typical set of parameters e — 600 MeV, 7=1 GeV and 
Oc = 0.1 GeV^, and we find weak dependence of the expo- 
nents on these parameters. That is, for the chain model-/, we 
find that the lowest energy chain state has higher energy than 
the bare state. In the critical case the energy increases less 
rapidly than the length of the chain, R, and is proportional to 
R^ for r>2. The average number of gluons grows weakly 
with R. The results are summarized in Figs. (|7]i, (jS} . 
In the case of model-// for r > 2, one easily finds, 

E[l=lir^l)yR[V,.H^Vm 
(N)" oc /;i/3 (72) 



I . , . , . , . , . 1 

20 40 60 80 100 

R [GeV '] 

FIG. 8. Expectation value of the the number of gluons in 
the chain as a function of QQ separation, R. A power law fit 
gives (A')^ = 0M5{R GeV)° ''23 and (N)" = 0.984(/? GeV)'"'*^ for 
model-I (solid line) and model-II (dashed line), respectively. 

while for r < 2 with the off-diagonal term dominating, 

E'o' = "lir^ l)YR[V,H-V,m- (73) 

Finally for the critical choice r = 2 numerical digitalization 
yields 

for the set of parameters, nig — 600 MeV, dc = 0. 1 GeV^, and 
7=1 GeV, Vc{°°) — Vc{0) — 1 GeV, and the results are shown 
in Fig. ^, (|9|. 

In model-// as /^-increases at some point the energy of the 
ground state chain increases less than the Coulomb potential. 
The chain state, however, has energy which is higher than 
that of the bare state, with the latter approaching a constant 
at large-/?. Thus in both models interactions among the chain 
increase the energy of the QQ pair as compared to the state 
with no gluons. 

IV. SUMMARY AND OUTLOOK 

We investigated microscopic origins of the gluon chain 
model. By analyzing the physical gauge interactions among 
constituent gluons, we fund a scenario for generating a chain. 
In this scenario a state with a number of gluons in the chain 
that is increasing with the separation between the QQ source 
emerges from interactions of dynamical gluons with the back- 
ground field. The background field is necessary in a phe- 
nomenological model of confinement if the latter is to orig- 
inate from condensation of chromomagnetic charges. These 
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FIG. 9. Ground state energy (solid line) of the chain Hamiltonian 
in model-//. A power law fit yields = 1.867(/?GeV)"-^^^GeV. 
The dashed line gives the energy of the bare state using for the string 
tension dc = 0.1 GeV^. 



interactions introduce off-diagonal elements into the effective 
Hamiltonian, which is one of the main differences between 
this and the chain model where the pair-production is absent. 
We have shown that the resulting ground state energy is con- 
vex ||32]| but the two models considered are still too simplistic 
to generate the linearly rising potential. While this deficiency 
can potentially be improved by considering more sophisti- 
cated models for the background field we found it difficult to 
reproduce the Zwanziger conjecture of "no-confinement with- 
out Coulomb confinement" [16 1 . We find the energy of the 
chain state to be higher then that of the bare one, defined as 
the expectation value of the Coulomb kernel in a state with no- 
backward reaction from the sources on the vacuum. It is possi- 
ble that a resolution of this problem requires renormalization 
for the single-gluon energies in the presence of the chain so 
that effectively e decreases with the number of gluons. 
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Appendix A: Background field model 



The correlation function G defined in Eq. ( 57 1 is computed 
using a classical distribution of sources of the background 
field. For example if these are monopole- antimonopole pairs 
the density p7v(c,c,) depends on the locations of the pairs. 
The expectation value of a function of is computed from 



(A?(x)A^'(y)) 



Z[A,A,] 
Z[l] 



(Al) 



where 
Z 



[^[Ac]] ^Zfdn" [ U^dCidCiP (c,c,) ^[A,] (A2) 



and the background field is given by 

N 

A^'(x) = n" ^[A„,(xx -Cx,/) - A„,(xx -Cx,;)] 
1=1 



(A3) 



where A„, is the abelian monopole field, and n" represents the 
(common) orientation of monopoles in the SU{Nc) algebra. 
For a uniform distribution of monopole-antimonopole pairs 
along the QQ-axis (z-axis) with the density given by 



P(c;,c,-) 



(A4) 



the background field is approximately constant along the QQ 



axis. In Eq. ( A4 1 p is the density of monopoles which is equal 
to the density of antimonopoles 



Px(jVg-l) 
R 

For the correlation function G{xj_,x) we then obtain 



G(xx,jc) = G(xx) =Px 



(A5) 



c/CxA„,(xx -Cx)Am(Cx 

^ yc/cxiiCxA„(cx)A,„(cx) 



(A6) 



The last term originates from the charge neutrality of the 
monopole-antimonopole distribution. If the core of the 
monopole field is smoothed out over a distance scale a = 
Oi^gho) then 



G(xx) --Pxlog 



|xxl 



(A7) 
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where the log/? dependence comes from cutting off the long 
range integral over the transverse plane. This is the standard 
expression for the correlation function of a pair of 2D vor- 
tices separated by a distance R. The 2D reduction originates 
from the assumption the monopoles are uniformly distributed, 
Eq. ( A5 I, along the direction of the QQ separation. 
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